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ABSTRACT 

We adapt an existing data compression algorithm, MOPED, to the extraction 
of median-likelihood star formation histories from medium-resolution galaxy spectra. 
By focusing on the high-pass components of galaxy spectra, we minimize potential 
uncertainties arising from the spectro-photometric calibration and intrinsic attenu- 
ation by dust. We validate our approach using model high-pass spectra of galaxies 
with different star formation histories covering the wavelength range 3650-8500 A at 
a resolving power of ~ 2000. We show that the method can recover the full star for- 
mation histories of these models, without prior knowledge of the metallicity, to within 
an accuracy that depends sensitively on signal-to-noise ratio. The investigation of the 
sensitivity of the flux at each wavelength to the mass fraction of stars of different ages 
allows us to identify new age-sensitive features in galaxy spectra. We also highlight a 
fundamental limitation in the recovery of the star formation histories of galaxies for 
which the optical signatures of intermediate-age stars are masked by those of younger 
and older stars. 

As an example of application, we use this method to derive average star formation 
histories from the highest-quality spectra of typical (in terms of their stellar mass), 
morphologically identified early- and late-type galaxies in the Early Data Release 
(EDR) of the Sloan Digital Sky Survey (SDSS). We find that, in agreement with the 
common expectation, early-type galaxies must have formed most of their stars over 
8 Gyr ago, although a small fraction of the total stellar mass of these galaxies may be 
accounted for by stars with ages down to 4 Gyr. In contrast, late-type galaxies appear 
to have formed stars at a roughly constant rate. We also investigate the constraints 
set by the high-pass signal in the stacked spectra of a magnitude-limited sample of 
20,623 SDSS-EDR galaxies on the global star formation history of the Universe and 
its distribution among galaxies in different mass ranges. We confirm that the stellar 
populations in the most massive galaxies today appear to have formed on average 
earlier than those in the least massive ones. Our results do not support the recent sug- 
gestion of a statistically significant peak in the star formation activity of the Universe 
at redshifts below unity, although such a peak is not ruled out. 

Key words: methods : data analysis - galaxies: statistics - galaxies: evolution - 
galaxies: stellar content 



1 INTRODUCTION 

Modern spectroscopic galaxy surveys suc h as the 2dF 
Galaxy Redhift Survey JColless et ail 1200 lT) and the Sloan 
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Digital Sky Survey (SPSS. Istoughton et al"ll2002l) are col- 
lecting hundreds of thousands of medium-resolution spec- 
tra of galaxies in the nearby un iverse. At higher red shift, 
the VIRMOS-VLT Peep Survey ihe Fevre et al.ll200ll) and 
the P eep Extragalactic Evolutionary Probe JPavis et all 
2003) are collecting the spectra of nearly 200 000 galaxies in 
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the younger universe. The high-quality spectra gathered by 
these surveys exhibit a myriad of stellar absorption features 
containing valuable information about the physical prop- 
erties of the galaxies, such as age, star formation history, 
metallicity and dust content. The extraction of this infor- 
mation from medium-resolution galaxy spectra is essential 
for improving current constraints on galaxy evolution. 

Stellar population synthesis, the modeling of the light 
emitted by specific populations of stars, is a natural ap- 
proach to interpreting observed galaxy spectra in terms of 
physical parameters. Until recently, the spectral resolution 
of standard population synthesis codes was typically 3-10 
times lower than that achieved in modern spectroscopic 
galaxy surveys (A/ A A ~ 2000 for the SDSS). Thus, for spec- 
tral analyses to be performed, observed high-quality galaxy 
spectra had to be degraded to the resolution of the mod- 
els, resulting in the loss of valuable high-frequency informa- 
tion. The advent of new models with higher spectral resolu- 
tion has opened the door to more refined spectral analyses 
jVazdekisHl999l : lBruzual fc Charlodl2003l) . For example, the 
ability to model in detail the 4000 A discontinuities and H<5a 
absorption indice s of m ore than 10 5 SDSS galaxies allowed 
iKauffmann et all ((2003) to derive useful constraints on the 
mass fraction of stars formed in recent bursts in these galax- 
ies. The next step is to identify the information provided by 
the entire medium-resolution spectrum on the star forma- 
tion history, metallicity and dust content of a galaxy. 

The extraction of physical parameters from large num- 
bers of medium-resolution galaxy spectra requires the de- 
velopment of dedicated techniques. Several algorithms have 
been designed to efficiently 'compress' galaxy spectra into 
a reduced number of parameters with minimum loss of in- 
formation. Such alg orithms include, fo r example, principal 
com ponent analysis ( Ronen et alJl999l) , Fisher matrix anal- 
vsis llTegmark et alJ 1997ft and information bottleneck (see 
lLahavi200ll for a review). Some of these algorithms have also 
been used to classify galaxy spectra, complementing pure 
classification methods such as trained neural networks. So 
far, however, these methods have not been used to extract 
physical parameters from medium-resolution galaxy spectra. 

In this paper, we use an original approach to extract 
star formation histories (and rough constraints on the metal- 
licities) from medium-resolution galaxy spectra. Our ap- 
proach relies on the existing algorithm MOPED 1 (Multi- 
ple Optimised Paramet er Estima t ion a nd Data compres- 
sion) developed bv lHeavens et ail <l200Cf) . which allows one 
to compress galaxy spectra into a reduced number of linear 
combinations connected to 'physical parameters' (for exam- 
ple, the mass fraction of stars in a given age range). MOPED 
has been used in the past to recover star formation his- 
tories, metallici ties and dust contents from low-resolutio n 
galaxy spectra dReichardt et al-lhoOll : Ipanter et alj [2003). 
The originality of our work consists in combining MOPED 
with a rece nt stellar population synt hesis model of medium 
resolution feruzual fc CharlotlEool hereafter BC03). We 
focus on the high-pass components of galaxy spectra, which 
allows us to minimize potential uncertainties arising from 
the spectro-photometric calibration and from the intrinsic 
attenuation by dust. This requires us first to remove any con- 

1 The MOPED algorithm has a patent. 



tamination of the high-frequency signal by nebular emission 
lines. Our approach also differs from previous applications 
of the MOPED algorithm in that we compute full likelihood 
distributions of physical parameters and do not rely on ac- 
celeration techniques to selectively explore the parameter 
space. This allows us to construct more realistic confidence 
intervals for the derived star formation histories than has 
been previously possible. 

We show that our method succeeds in recovering the 
star formation histories from a variety of model galaxy spec- 
tra, to an accuracy that increases with signal-to-noise ratio. 
We investigate the spectral features that are most sensi- 
tive to stars in various age ranges, excluding those regions 
of the optical spectrum that are known to depend strongly 
on changes in metal abundance ratios. We also point out 
a fundamental limitation in the spectral interpretation of 
galaxies for which the optical signatures of intermediate-age 
stars are masked by those of younger and older stars. To 
illustrate the usefulness of our approach, we use it to inter- 
pret the spectra of galaxies with high signal-to-noise ratios 
in the SDSS early data release (EDR). We find that, in gen- 
eral, galaxies with high concentrations typical of early-type 
galaxies appear to have formed the majority of their stars 
prior to z = 1, while late- type, low-concentration galax- 
ies appear to have formed stars at a roughly constant rate 
at all times. From the analysis of the stacked spectra of 
a magnitude-limited sample of 20,623 SDSS-EDR galaxies, 
we can also investigate the global cosmic star formation 
history of the Universe. Our results confirm previous evi- 
dence that the stars that end up in massive galaxies today 
formed on average earlier than those that end up in low-mass 
galaxies ()Heavejrs_et_alJ |2004|; see also ICowie et al.l Il996l : 
iKodama et alJl2004l : IKauffmann et al.ll2004l) ~ Our analysis 
does not reveal any statistically significant peak in the in- 
tegrated star formation activity of the Universe at redshifts 
less than z ~ 1. 

The paper is organized as follows. In Section |^1 we de- 
scribe our approach and the population synthesis models we 
rely on. We also recall the main features of the MOPED al- 
gorithm. In Section |3J we use model simulations to evaluate 
the extent to which star formation histories can be recovered 
from medium-resolution galaxy spectra. We also investigate 
the spectral features that are most revealing of stars in var- 
ious age ranges. In Sectional we first derive star formation 
histories from individual, high-quality spectra of high- and 
low-concentration galaxies in the SDSS EDR. We also in- 
vestigate the typical star formation histories of high- and 
low-concentration galaxies, based on the analysis of stacked 
spectra. We complement this with a study of the overall 
star formation history of the Universe using a sample of 
20,623 galaxies from the SDSS EDR. Our conclusions are 
summarized in Section |S] Throughout this paper, we adopt 
a Friedmann-Robertson- Walker cosmology with Q — 0.3, 
A = 0.7 and H = 70kms~ 1 Mpc" 1 . 



2 THE METHOD 

We summarize here our approach for extracting physical pa- 
rameters from medium-resolution galaxy spectra. We first 
briefly review the main features of the compression algo- 
rithm MOPED and how we adapt it to interpret the high- 
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pass signal in galaxy spectra. We describe the physical pa- 
rameters of interest to us and how we estimate their likeli- 
hood. 



plication of this algorithm to the interpretation of medium- 
resolution galaxy spectra. 



2.1 The MOPED algorithm 

MOPED is a linear compression algorithm, which effi- 
ciently reduces a large array of data points into an ar- 
bitrarily small number of key parameters. Such an algo- 
rithm is ideally suited to the extraction of physical param- 
eters from medium-resolution galaxy spectra . We refer the 
reader to the o riginal paper of lHeavens et all ( 2000, sec also 
iTeemark et"ai]|l997l) for a more detailed description than 
given in the following summary. 

• In MOPED, the compression of a galaxy spectrum F\ 
(defined as a vector of fluxes per unit wavelength) into a 
small number np of physical parameters - for example, the 
mass fraction of stars in a given age range - is represented 
by a vector "y* of np elements yj, defined by yj = T> j • F\. 

• The np weighting vectors hj are specific to the spec- 
trum Fx- They quantify in an optimal way the sensitivity of 
the flux at each wavelength to variations in the param- 
eter, taking into account the noise associated to F\. The 
construction of the weighting vectors b j requires a popula- 
tion synthesis model to compute the derivatives of the flux 
at each wavelength with respect to each physical parameter. 
In practice, for the sake of efficiency, these derivatives are 
estimated at a single point in the parameter space called 
the 'fiducial model' (see Section [2.4. II below). The resulting 
vectors hj contain information about the spectral features 
most sensitive to the parameter j, given the noise pattern. 

• Once the weighting vectors bj are determined, the 
probability distributions of the np parameters given the 
spectrum Fx may be evaluated by computing the likelihood 
of each (compressed) model ^ nlod in the parameter space, 
according to the norm of the vector y* — y*mod- 

The MOPED compression would not cause any loss of 
information about the np parameters if these had Gaus- 
sian likelihood distributions, and if the fiducial model cor- 
responded for each F\ to the best fit model that would be 
obtained without compression. In most applications how- 
ever, these criteria are not met (Section 12. 4. 21 . It is also 
worth noting that, because the weighting vectors hj are 
obtained recursively in a process similar to a Gram-Schmidt 
orthonormalization of independent vectors, the values of the 
weights may depend on the order in which the parameters 
are considered. This is especially true when some spectral 
features present similar sensitivities to several parameters 
(see also Section f2.3.2ll . We have found that, in general, our 
results are relatively robust to changes in the order in which 
the weighting vectors b j are e valuated. We therefore adopt 
here the natural or der used bv lReichardt et all i200lf) and 
IPanter et alJ (|2003), which consists in evaluating the sen- 
sitivity of the spectrum to the oldest stars first, and then 
proceed with younger stars. 

As mentioned in Section the MOPED algorithm has 
already been used to extract star formation histories from 
low-resolution spectra of galaxies in the Kennicutt (1992) at- 
las an d the SDSS survey ( Reichard t et al ■l200lHPanter et alJ 
2003). In the remainder of this section, we describe our ap- 



2.2 Medium-resolution galaxy spectra 

To interpret medium-resolution galaxy spectra with 
MOPED, we rely on the recent population synt h esis c ode 
of BC03. This is based on the iLe Borgne et al] (12003ft li- 
brary of observed stellar spectra and has a resolution of 3 A 
across the whole wavelength range from 3200 A to 9500 A 
(corresponding to a median resolving power A/AA w 2000 
and a nominal velocity dispersion of 70 km s" 1 ), for a wide 
range of metallicities (—2.0 < [Fe/H] < +0.50). Since we are 
interested here in the interpretation of SDSS galaxy spectra 
(see Sectional, we restrict our analysis to the wavelength 
range from 3650 to 8500 A, where the noise is lowest in the 
rest-frame SDSS spectra. To speed up the analysis, we in- 
crease the wavelength sampling of the BC03 model spectra 
from the original 1-A sampling to a 3-A sampling. Further- 
more, in all applications below, we always use model spectra 
broadened to a velocity dispersion consistent with that of the 
galaxy under study. 

The observed stellar spectra incorporated in the BC03 
models come from nearby stars with solar abundance ra- 
tios. However, stars in external galaxies are known to ex- 
hibit changes in heavy-element abundance ratios. In partic- 
ular, the ratio of a elements to iron appears to be enhanced 
in massive early-type galaxies relative to the solar com po- 
sition (e.g.. lWorthev et al]ll99l lEisenstein et aljl2003l) . To 
increase the applicability of our models to the interpretation 
of observed galaxy spectra, we choose to mask out those re- 
gions of the optical spectrum that are known to strongly 
depend on changes in metal abundance ratios. 

BC03 used a large library of model star formation his- 
tories to fit strong atomic and molecular features in the 
highest-quality spectra of the SDSS EDR. They consid- 
ered the 25 commonly measured spectral features of the 
'Lick system' that were defined and calibrated in the spec- 
tra of 460 Galactic stars in the wavelength range from 
4000 A to 6400 A jWorthevI Il994l: IWorthev fc Ottavianil 
Il997t iTrager et alJ ll99^) ~As expected, the Lick features 
that are the least well fitted by the BC03 models in the 
SDSS spectra are those Frhomas et al] i2003T) found to de- 
pend most strongly on changes in a/Fe. We therefore ex- 
clude from our analysis the spectral b andpasses corre s pond - 
ing to the Lick features found by iThomas et al] (120031) 
to be most sensitive to changes in a/Fe. These are Mgi, 
Mg 2 , Fe4383, Fe5335 and Fe5406. The wide TiOi and 
Ti02 features are also not always well reproduced by the 
BC03 models in the highest-quality SDSS spectra. How- 
ever, we keep them in the present analysis because they 
are thought to depend less strongly on a/Fe (iThomas et al] 
2003). We also do not exclude the spectral bandpasses defin- 
ing other Lick indices, which sometimes overlap with the 
bandpasses of Mgi, Mg 2 , Fe4383, Fe5335 and Fe5406. In 
summary, the spectral windows we exclude from our analysis 
are 4363-4424 A, 4445-4458 A, 4897-4958 A, 5071-5135 A, 
5157-5199 A, 5318-5367 A and 5377-5425 A, corresponding 
to a total spectral width of 338 A. 

Our main motivation in the present study is to identify 
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the constraints set by the high-pass signal 2 in galaxy spec- 
tra on the past history of star formation. We choose not to 
include the low-pass components of galaxy spectra for sev- 
eral reasons. Firstly, the low-pass signal is strongly affected 
by dust attenuation and potential spectro-photometric in- 
accuracies, which introduce unwanted uncertainties in the 
derived star formation histories. Secondly, population syn- 
thesis models with medium spectral resolution have only 
recently become available, and it is important to establish 
the extent to which the high-pass signal helps us constrain 
the physical parameters of galaxies. Finally, the MOPED 
algorithm has already been used to interpret low-resolution 
galaxy spectra in terms of star formation histories (see Sec- 
tion iiuj. 

We compute the high-pass component F^ F of a galaxy 
spectrum Fx (either modelled or observed) by first smooth- 
ing Fx with a top-hat kernel W(X) of smoothing diameter 
l s = 500 A. This yields the low-pass spectral component 



A+is/2 



dX' W(X — A') Fy 



(1) 



A-i s /2 

The high-pass component may then be obtained as either 
Fa 



F A HP — Fx — F A LP . 



(2) 



(3) 



Equation @, which was used by iBaldrv et all (2002), pre- 
serves the equivalent widths of individual spectral features. 
There is no strong a priori motivation for using equations 
or to compute the high-pass components of galaxy spec- 
tra. Equation which would be preferred if we achieved 
true frequency filtering, appears to perform better in the re- 
covery of star formation histories in model simulations. We 
therefore use it throughout this paper to compute high-pass 
galaxy spectra (we have checked that using equation|!2]would 
not affect our main conclusions). Since the noise associated 
to the low-pass component is negligible relative to that asso- 
ciated to the high-pass component, the signal-to-noise ratio 
per pixel in the high-pass spectrum is virtually the same as 
that in the original spectrum. 

The high-pass galaxy spectra obtained using equa- 
tion © must be normalized for the purpose of comparisons 
between models and observations. Ideally, the normalization 
should be achieved in a band which traces the total stellar 
mass of the galaxy. In practice, we are limited by the wave- 
length coverage and redshifts of the SDSS spectra we plan 
to analyze. We have found that the rest-frame 6000-7300 
A band provides a good compromise and therefore normal- 
ize all spectra below in this wavelength rang e. All results 
presented here are based on a IChabrierl J2003T) initial mass 
function with lower and upper cutoffs 0.1 and 100 M . 



2 In this paper, we use the term 'high pass' to designate that part 
of a galaxy spectrum containing high-frequency information. We 
prefer the term 'high pass' over 'high frequency' because we do 
not strictly perform true frequency filtering. 



2.3 Physical parameters 

We now describe the physical parameters we wish to recover 
from medium-resolution galaxy spectra. We are primarily 
interested in the past history of star formation. We also de- 
rive constraints on the stellar metallicity and mention the 
potential of our approach for constraining the attenuation 
of starlight by dust in galaxies. 



2. 3. 1 Star formation history 

We characterize the star formation history of a galaxy by 
constraining the evolution of the star formation rate in 6 
consecutive age bins. Ideally, these age bins should corre- 
spond to significant changes in the spectral characteristics 
of the stars, such as those arising through the onset of red su- 
pergiants, bright asymptotic giant stars, late-B and A stars 
with strong H-Balmer absorption lines, and red giant branch 
stars of different temperatures. 3 After some experimenta- 
tion, we settled on the following age bins to characterize 
the star formation history of a galaxy over the period from 
the age of the Universe to the present time: [13.5-8], [8-4], 
[4-1.5], [1.5-0.5], [0.5-0.1] and [0.1-0] Gyr. Within each bin, 
star formation is assumed to proceed at a constant rate. The 
relative 'amplitude' of a bin then determines the fraction of 
the total stellar mass of a galaxy that formed at the epoch 
corresponding to that bin. 

The sampling of galaxy star formation histories in the 
MOPED algorithm is defined by the possible combinations 
of star formation amplitudes in the above 6 age bins. We 
have chosen 6 levels of star formation amplitude in each 
age bin, such that the fraction of the total stellar mass of 
a galaxy formed in any age bin be sampled in an optimal 
way among all possible star formation histories. For each of 
the 3 oldest age bins, we require that this fraction be lin- 
early well sampled (i.e. with approximately equal numbers 
of models) by 6 values between and 1. For each of the 3 
youngest age bins, we require that it be logarithmically well 
sampled by 6 values between 10~ 3 and 1 (we do not dis- 
tinguish fractions of 10~ 3 or less from zero star formation). 
To select the absolute star formation amplitudes providing 
the closest approximation to such a sampling in stellar mass 
fractio n, we have used the non -negative least squares rou- 
tine of lLawson fc Hanson! <I1974T) . We find that 6 amplitudes 
per bin are sufficient to sample the parameter space in an 
adequate way. Adopting a computationally more expensive 
sampling of 8 amplitudes per bin does not significantly im- 
prove the accuracy to which the method can recover star 
formation histories from medium-resolution galaxy spectra 
(Section 123. 



3 iReichardt et alj <200ll) and lPanter et alJ <2003ft use 8 age bins 
to describe the star formation history of a galaxy in their analysis 
of low-resolution spectra using MOPED. While the computation- 
ally more demanding interpretation of medium-resolution spectra 
prevents us from using such a large number of bins, we find that 
adopting more than 6 bins does not statistically improve the con- 
straints we derive on the star formation history. This is because 
the spectral signatures of stars of too similar ages cannot be easily 
distinguished, even at medium resolution. 
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2.3.2 Metallicity 

For the purpose of the present paper, we assume that stars 
of all ages in a given galaxy have the same metallicity. We 
interpret this as the 'effective' stellar metallicity (i.e., the 
metallicity of the stars dominating the light). This choice is 
motivated primarily by the fact that our models have solar 
heavy-element abundance ratios, while the abundance ratios 
of he avy elements are known to vary in extern al galaxies 
(e.g., IWorthev et~al]ll992t lEisenstein et al.ll2003ft . As men- 
tioned in Section fe.^l we exclude from our analysis the spec- 
tral features with the strongest identified dependence on 
changes in the a/Fe abundance ratio. While we believe that 
this is sufficient for obtaining rough constraints on effective 
metallicities, we are less confident in the reliability of the 
models for recovering accurate metal enrichment histories 
from observed galaxy spectra. This caution is further jus- 
tified by the fact that the signatures of metals are often 
similar in the spectra of stars with different ages. Allowing 
different metallicities for different age bins would therefore 
introduce an unwanted dependence of the results on the or- 
der in which the MOPED weighting vectors are computed 
(see Section \'2. Q . We note that the star formation histories 
derived in Sections [3] and [I] below are marginalized over ef- 
fective metallicity (Section 12. 4. 21 . 

For simplicity, we consider here only 3 choices for the 
effective stellar metallicity of a galaxy: Z — OAZq, Zq and 
1.5Z@. These values are representative of the range in stellar 
metallicities inferred from the analysis of selected features 
with negligible dep endence on a/Fe in t he spectra of ~ 2 x 
10 s SDSS galaxies JOallazzi et al.ll2004) . 



2.3.3 Attenuation by dust 

The analysis presented in this paper, which relies on the in- 
terpretation of high-pass optical galaxy spectra, is largely 
exempt from the uncertainties arising from attenuation of 
starlight by dust in galaxies (Section 12.21 . The potential 
weak effects of dust on the relative strengths of stellar ab- 
sorption features are expected to affect only negligibly the 
star formation histories derived from high-pass galaxy spec- 
tra (see Section f3.2H . 

It is worth mentioning that our approach also provides 
a means of constraining attenuation by dust. The star for- 
mation history (and effective metallicity) derived from the 
high-pass spectrum of a galaxy can be exploited to predict 
the low-pass spectrum of that galaxy using dust-free pop- 
ulation synthesis models. By comparing this low-pass spec- 
trum with the observed one, it is possible to derive useful 
constraints on the attenua tion by dust in the ga laxy (for ex- 
ample, using the model of lCharlot fc Falll fcoOOl. Unlike the 
analysis of the high-pass spectrum, this requires an accurate 
spectro-photometric calibration of the observed spectrum. 



2.4 Statistical estimates of physical parameters 

One of the main advantages of the MOPED algorithm over 
other techniques to interpret galaxy spectra is that it allows 
the derivation of statistical estimates of physical parameters. 
We outline below the way in which this is achieved in our 
models. 



2.4-1 Choice of the fiducial model 

In the MOPED algorithm, the fiducial model is the point of 
the parameter space where the weighting vectors describing 
the sensitivity of the spectrum to the various physical pa- 
rameters are evaluated (Section 12. H . The vectors are con- 
structed by computing the derivatives of the flux at each 
wavelength with respect to each parameter and accounting 
for the noise pattern. 

The requirement in choosing the fiducial model is that 
its spectrum should contain as many spectral features as 
possible to constrain the different parameters of interest. To 
constrain the star formation history, the spectrum should 
therefore have clear signatures of stars in all age ranges. In 
the present study, the choice of the fiducial model is simpli- 
fied by the fact that we are interested only in the high-pass 
components of galaxy spectra, as obtained by subtracting 
the low-pass components from the total spectra (equation^ ■ 
Also, in our models, galaxy spectra are simple linear com- 
binations of the spectra of stellar populations of identical 
metallicity in 6 age bins (Section l2.3.H . Finally, we neglect 
any potential weak, non-linear effect of dust on high-pass 
galaxy spectra. All this implies that, at fixed metallicity 
and stellar velocity dispersion, the derivatives of the high- 
pass component with respect to physical parameters do not 
depend on the specific model (as characterized by the ampli- 
tudes of the stellar populations in the 6 age bins) at which 
they are evaluated, so long as this model includes contribu- 
tions from all age bins. 4 We compute such a set of high-pass 
derivatives for each of the 3 metallicities considered here, 
Z = 0.4^0, Zq and 1.5Zq, and for several stellar velocity 
dispersions between 70 and 300 km s _1 (Sections 12.21 and 



We take the noise spectrum entering the computation 
of the MOPED weighting vectors at the fiducial model to be 
that of the observed galaxy under study. The disadvantage 
is that this requires us to evaluate the weighting vectors for 
every galaxy in the observed sample. 5 The reason for this 
refinment is to limit the loss of information about the phys- 
ical parameters ar ising from the compre ssion of the spectra. 
As emphasized bv lHeavens et alJ feOOOf), the MOPED com- 
pression is optimal when the properties of the fiducial model 
are as close as possible to those of the observed spectrum. 

2.4-2 Probability density functions 

The likelihood that any model in the parameter space repre- 
sents a given observed spectrum may be computed by com- 
paring the compressed model and observed spectra as de- 
scribed in Section [2. II If the likelihood distributions of the 
physical parameters were well approximated by multivariate 
Gaussian functions over the entire parameter space, sam- 
pling the neighbourhood of the most likely model would be 

4 It is interesting to note that this would not be the case if we had 
chosen equation instead of equation J3J to compute high-pass 
galaxy spectra. In this case, the low-pass component would enter 
the expression of the derivatives of the high-pass component with 
r espect to the phy s ical p aram eters. 

5 iReichardt et alJ <200ll) and lPanter et alJ <2003ft adopt a fixed, 
'average' noise spectrum to interpret low-resolution galaxy spec- 
tra with MOPED. 
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sufficient to estimate the confidence intervals for these pa- 
rameters. In our analysis, we have found that the likelihood 
distributions of the fractions of stars formed in different age 
bins is not necessarily well approximated by Gaussian func- 
tions. Thus, the parameter space must be explored more 
thoroughly to derive constraints on the star formation his- 
tories of galaxies. 

Sampling the entire parameter space to evaluate the 
likelihood of every model given an observed spectrum is 
computationally demanding. For this reason, previous stud- 
ies appealed to standard techniques such as 'Monte-Carlo 
Markov Chains' to selectively explore the parameter space 
and const ruct the likelihood distrib utions of physical pa- 
rameters iPanter et al.ll20o3 : see also lVerde et alJl2003fl . In 
practice, such techniques rely on the sampling of the pa- 
rameter space by a fixed, limited number of points. As a 
consequence, it is unclear whether the parameter space is 
always sampled in an adequate way to estimate the true 
probability density functions of physical parameters. 

We choose here to favor the accuracy of error estimates 
over the reduction of computational cost and always explore 
the entire parameter space to compute the likelihood distri- 
butions of the different physical parameters. We first com- 
pute, for each metallicity, the Bayesian likelihood distribu- 
tion of the fraction of stars in each age bin. This is evalu- 
ated by weighting each model by its probability and then 
summing the probabilities in bins of fraction of stars (see 
appendix A of [Kauffmann et al.ll20oj) . We then marginalize 
the likelihood distribution of the fraction of stars in each age 
bin over metallicity, by summing the likelihood distributions 
over the 3 metallicities Z — OAZq, Zq and 1.5Zq. In this 
procedure, we assume that the 3 metallicities are equally 
probable a priori. We refer to the result of this sum as the 
probability density function of the fraction of stars in that 
age bin. By construction, this function includes the errors 
associated to the metallicity determination. We also obtain 
a rough estimate of the most likely metallicity of the galaxy 
by marginalizing the likelihood distributions obtained for 
each metallicity over star formation histories. 



2.4.3 Priors 

In our models, the prior distributions of the parameters are 
flat, reflecting a point of view of 'maximum ignorance'. The 
fraction of stellar mass formed is uniformly distributed in 
each of the oldest 3 bins (ages ^ 1.5 Gyr) and logarithmically 
distributed in each of the youngest 3 bins (ages < 1.5 Gyr; 
see Section fe.B.H . In addition, the models are uniformly dis- 
tributed in metallicity between Z = OAZq, Zq and 1.5Zq 
(Sections 12.3.21 and 12.4.21 . We find that our results are ro- 
bust to changes in these prior distributions. For example, 
changing from a uniform to a logarithmic sampling of the 
stellar mass fractions (and conversely) has little effect on 
the derived star formation histories. Such an invariance is 
expected when the problem is well constrained. 



3 ASSESSMENT OF EFFICIENCY USING 
MODEL GALAXY SPECTRA 

3.1 Outline 

This Section is the core of the paper. Here, we use model 
simulations to assess the accuracy to which our method 
is able to recover star formation histories from medium- 
resolution galaxy spectra. For illustration, we parametrize 
the star formation history as a simple exponential law 
SFR(t) oc exp (— t/r), where t is the age of the galaxy and r 
the time scale of star formation. We consider 3 models with 
t = 2, 7 Gyr and 00 (i.e. constant star formation rate), all 
with the same present galaxy age t — 13.5 Gyr (correpond- 
ing to the age of the Universe for the adopted cosmology). 
These simple models are often taken to represent the star 
formation hist ories of, respectively, E/SO, S b/Sc and Sm/Im 
galaxies (e.g., iKauifmann fc Charlo1ill998i) . For simplicity, 
we assume that all 3 models have solar metallicity and a 
stellar velocity dispersion cry — 100 km s _1 (we discuss be- 
low the effects of changing these parameters). We consider 
two values of the median signal-to-noise ratio per pixel in 
the spectra, S/N = 10 and S/N = 30. To be representative, 
the results presented here are always averaged over 50 in- 
dependent realizations of the noise (increasing this number 
from 50 to 100 has a negligible effect on the results). 

In Section 13.21 below, we evaluate how accurately our 
method can recover the star formation histories of these 
models, based on the analysis of high-pass spectra. Then, 
in Section ^. 31 we point out a fundamental limitation in the 
recovery of the star formation histories of galaxies for which 
the optical signatures of intermediate-age stars are masked 
by those of younger and older stars. Finally, in Section f3. 41 
we highlight the spectral features that are the most sensitive 
to the presence of stars in various age ranges. 



3.2 Recovery of star formation histories 

We first consider the case of a model galaxy with a star for- 
mation time scale r = 7 Gyr. This model is particularly in- 
teresting, because it corresponds to the expected time scale 
of star formati on of the Universe when averaged over large 
volum es (e.g., iGlazebrook et all |2003; Br inchmann et all 
2004). We adopt for the moment S/N = 30 and take the 
wavelength dependence of the noise to be that observed for 
a typical late-type galaxy in the SDSS (Section 14. H . 

Fig- shows the star formation history recovered from 
the high-pass spectrum of this galaxy using the method pre- 
sented in Section [5] Each panel shows the average probabil- 
ity density function of the mass fraction of stars formed in 
a given bin of look-back time. This was obtained by sum- 
ming the 50 probability density functions obtained for dif- 
ferent realizations of the noise (see above). The results are 
shown in a more compact way in Fig. where we sum- 
marize the probability distribution in each age bin by the 
median value (thick line) and 16-84 (box) and 2.5-97.5 (er- 
ror bars) percentile ranges. We stress that these ranges do 
not reflect the scatter in the median mass fractions recov- 
ered from 50 different realizations of the noise. Instead, for 
each age bin, the boxes and error bars in Fig. indicate 
the 68 and 95 percent confidence intervals measured from 
the corresponding average probability density function of 
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Figure 1. Test recovery of the star formation history from the 
high-pass spectrum of a 13.5 Gyr old stellar population with a 
star-formation time scale of 7 Gyr and the solar metallicity. In 
each panel, the histogram shows the average probability density 
function of the mass fraction of stars formed in a given bin of 
look-back time (as indicated). This was obtained by summing the 
probability density functions obtained for 50 different realizations 
of the noise. The dotted line shows the true fraction of stellar 
mass formed in the original model, while the dashed line shows the 
most likely model recovered by the analysis. The median signal- 
to- noise per pixel in the input spectrum is S/N = 30, and the 
stellar velocity dispersion cry = 100 km s _1 . 



Fig. □ These ranges would correspond to the ±1<t and ±2<r 
confidence intervals in the case of Gaussian distributions. 
For reference, the dotted lines in Fig. Q and the triangles in 
Fig. [5^ show the true fractions of stellar mass formed at the 
different look-back times in the input model. We recall that 
the probability distributions in Figs0and|2Jt include the un- 
certainties associated to the assumed unkown metallicity of 
the stars (Section 12. 4. 2H . The solar metallicity of the model 
is recovered with 70 percent confidence. 

The results of Figs and indicate that, overall, the 
method recovers remarkably well the star formation history 
of the test galaxy. The true mass fraction of stars formed at 
different look-back times is always recovered within the 16- 
84 percent (±lcr) confidence range predicted by the analysis. 
Some discrepancies exist between the recovered and true star 
formation histories. In particular, the mass fraction of stars 
with ages around 1 Gyr appears not to be well constrained. 
Its median value is underestimated, while the median mass 
fractions of stars in the adjacent age bins are overestimated. 
In Section 13.31 below, we show that this mismatch arises 
from a general limitation of the spectral interpretation of 
galaxies for which the optical signatures of intermediate-age 
stars (0.5-4 Gyr) are masked by those of younger and older 
stars. 

The accuracy to which the star formation history can be 
recovered from high-pass galaxy spectra depends sensitively 
on the signal-to-noise ratio. Fig. |5Jd shows the same results 
as in Fig. Ha,, but assuming a median signal-to-noise ratio 
per pixel of 10 instead of 30 in the input spectrum. In this 
case, the agreement between the predicted and the true star 
formation histories is significantly worsened. For example, 
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Figure 2. (a) Same as Fig. but shown in a more compact 
way. In each bin of stellar age, the thick line shows the median 
fraction of stars predicted to have formed, based on the analysis 
of the model high-pass spectrum, while the vertical size of the 
box and the error bars show the 16-84 and 2.5-97.5 percentile 
ranges in this fraction, respectively (these ranges indicate the 68 
and 95 percent confidence intervals measured from the average 
probability density functions of Fig. they would correspond to 
the ±ltr and ±2cr confidence intervals in the case of Gaussian 
distributions). The triangles indicate the true fractions of stellar 
mass formed in the model, (b) Same as (a), but adopting a median 
signal-to-noise per pixel of 10 instead of 30 in the input spectrum. 



the method cannot recover the differences in the relative 
fractions of stars in the oldest 3 bins. For S/N = 10, the 
results also become sensitive to the order chosen to construct 
the weighting vectors in the MOPED algorithm (Section|2Jl. 

The effect of increasing the velocity dispersion at fixed 
S/N is expected to be similar to that of decreasing S/N at 
fixed velocity dispersion. In practice, the amplitude of the 
effect depends on the type of stars dominating the light. For 
example, changing the velocity dispersion from our adopted 
100 km s _1 to 300 km s _1 has only a weak effect on the re- 
sults presented in Figs|5K and b. 

We have also checked that attenuation by dust affects 
negligibly the results presented in Fig.|21 We used the sim ple 
but physically motivated model of lCharlot fc Fall i2000f) to 
include the effects of attenuation by dust on the spectrum 
of the t = 7 Gyr model galaxy. We adopted a typical V- 
band effective absorptio n optical depth fy = 1 appropriate 
for SDSS galaxies jBrinchmann et al-lEooi) . for which we 
assumed standard fractions of 70 and 30 percent to arise 
from dust in giant molecular clouds an d in the diffuse in- 
terstellar medium dCharlot fc Fallll2000l) . For S/N=30, the 
median star formation history recovered from the high-pass 
spectrum of this model does not differ in any systematic way 
from that shown in Fig. [5^, the discrepancies between the 
two cases being less than 4 percent in all age bins. 

Figs |3] and 01 show further examples of the recovery of 
star formation histories from high-pass model galaxy spec- 
tra. Fig. [3] shows the results obtained for a galaxy with a 
constant star formation rate, while Fig. 0] shows those ob- 
tained for a galaxy with a star formation time scale of 2 Gyr. 
In the latter case, we take the wavelength dependence of the 
noise to be that observed for a typical early-type galaxy in 
the SDSS (Section U|TJ|. For S/N = 30, the method recovers 
quite well the star formation histories of both model galax- 
ies (Figs |2K and 0K). The true mass fractions of stars in 
all age bins are recovered within the 16-84 percent (±lcr) 
confidence ranges predicted by the analysis. We note that, 
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Figure 3. Same as Fig. [3] but for a galaxy with constant star 
formation rate. 
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Figure 4. Same as Fig. [5] but for a galaxy with a star-formation 
time scale of 2 Gyr. 



for constant star formation rate (Fig. EJi), the fraction of 
stars with ages around 1 Gyr is somewhat better constrained 
than in the case of the 7 Gyr time-scale model (Fig. ■ As 
before, reducing the signal-to-noise ratio from 30 to 10 de- 
grades the quality of the results, especially for the relative 
fractions of stars in the oldest age bins. The solar metallicity 
of the models is recovered with 70 percent confidence. 

It is worth pausing here to emphasize the significance 
of the results presented in Figs These results are the 
first to quantify the usefulness of the high-pass components 
of medium-resolution galaxy spectra for constraining star 
formation histories. We have shown that, for spectra with 
high-enough signal-to- noise ratio (S/N ~ 30 per pixel), 
the method outlined in Section [5] succeeds in recovering 
with reasonable accuracy the star formation histories (and 
metallicities) of galaxies with a wide range of time scales 
of star formation . These results greatly extend the work of 
Kauf fmann et all (|^003), who derived valuable but limited 
constraints on the recent star formation histories of galaxies, 
based purely on the strengths of the 4000 A discontinuity 
and H5a stellar absorption line. We find that the analysis 
of the high-pass spectrum over the whole wavelength range 
from 3650 to 8500 A provides useful constraints on the en- 
tire star formation history of a galaxy, along with some con- 
straints on the stellar metallicity. 

A main advantage of our method over other approaches 
to extract star formation histories from galaxy spectra is 



that it minimizes uncertainties linked to attenuation by dust 
and the spectro-photometric calibration (Sections 12.21 and 
12.3.31 . These uncertainties, which often must be included in 
a rudimentary way, have impaired most previous interpre- 
tation of galax y spectra using low-resolution models ( e.g., 
most recently. iReichardt et al]l200lt|Panter et al.ll2003h . In 
this context, the results presented in Figs should be 
taken as illustrative of the new type of constraints that can 
be set on the star formation histories of galaxies, thanks to 
the advent of medium-resolution population synthesis mod- 
els. It is worth re-emphasizing that, in the examples shown 
above, the simulated galaxy spectra have the same solar 
heavy-element abundance ratios as the models used to in- 
terpret them, while the abundance ratios of heavy elements 
are known to vary in extern al galaxies (e.g.. IWorthev et alJ 
Il992t lEisenstein et al1l2003l) . As mentioned in Section |2~21 
we exclude from our analysis the spectral features with 
the strongest identified dependence on changes in the a/Fe 
abundance ratio. This should limit the uncertainties in the 
star formation histories and metallicities derived from ob- 
served, high-pass galaxy spectra (Section 0J. 



3.3 The case of intermediate-age stars 

In the previous section, we noted the difficulty in recover- 
ing the mass fraction of stars with ages around 1 Gyr from 
the high-pass spectrum of a 13.5 Gyr-old model galaxy with 
a star-formation time scale of 7 Gyr (Fig. This also 
affected the constraints derived on the mass fractions of 
stars in the adjacent age bins. We show here that this dif- 
ficulty arises from a fundamental limitation of the spectral 
interpretation of galaxies for which the optical signatures of 
intermediate-age stars (0.5-4 Gyr) are masked by those of 
younger and older stars. 

Fig. shows the relative contributions by stars in dif- 
ferent age bins to the integrated spectra (including both the 
high-pass and low-pass components) of the 3 model galaxies 
studied in Section l3~^l with star-formation time scales r = 2, 
7 Gyr and co (constant star formation rate). For r = 2 Gyr, 
only stars older than 4 Gyr contribute significantly to the 
integrated light at any wavelength. If the signal-to-noise ra- 
tio is high enough, the signatures of stars with ages 4-8 
and 8-13.5 Gyr can be discerned in the integrated high-pass 
spectrum (Fig.|5Jl. For constant star formation rate, stars of 
any age always contribute to a major fraction of the inte- 
grated light in at least some wavelength range. Young stars 
dominate at short wavelengths, while intermediate-age and 
old stars dominate at long wavelengths. This makes the con- 
tributions by stars of all ages identifiable in the integrated 
high-pass spectrum of the galaxy (Fig. 

The case of the t — 7 Gyr model is different. The few 
young stars in this model dominate the integrated light at 
short wavelengths. However, the large number of old stars 
dominate so much the emission at long wavelengths, that 
intermediate-age stars (0.5-4 Gyr) never contribute more 
than a minor fraction of the integrated light at any wave- 
length. This is why the mass fraction of intermediate-age 
stars is difficult to constrain accurately from the integrated 
high-pass spectrum of the galaxy, even for very high signal- 
to- noise ratio (Fig. 0. This limitation does not affect the 
interpretation of only high-pass spectra. As Fig. shows, 
it is a fundamental limitation of the spectral interpretation 
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of galaxies with smoothly declining star formation rates, for 
which the signatures of intermediate-age stars are masked 
by those of younger and older stars (an analogous limitation 
can arise when the spectrum is dominated by large amounts 
of recent star formation; see Section 14,21 and Fig. 1111 be- 
low). We note that r » 7 Gyr is thought to be the volume- 
averaged time scale of star formation of the Un iverse (e.g., 
iGlazebrook et al ■l2003l : lBrinchmann et all2004Tl . Hence, our 
results imply that it may be difficult to constrain with pre- 
cision the star formation history of the Universe at redshift 
z ~ 0.5, based purely on the spectra of 2 ~ 0.1 galaxies 
(Section EH. 



3.4 New insight into spectral indicators 

It is instructive to investigate the spectral features selected 
by the MOPED algorithm to constrain the star formation 
histories of galaxies from high-pass spectra. Given a fixed 
noise pattern, the weighting vectors constructed as outlined 
in Section \'2. II quantify in an optimal way the sensitivity of 
the flux at each wavelength to the mass fractions of stars in 
various age ranges. Thus, these vectors provide a means of 
identifying optimal spectral indicators of the star formation 
histories of galaxies. 

Fig. [U shows the components of the weighting vectors 
quantifying the sensitivity of the flux at each wavelength 
to the mass fractions of stars in the 6 age bins used in our 
analysis (red histogram) . In this example, we have adopted a 
stellar velocity dispersion ov = 100 km s~\ solar metallicity 
and a fixed signal-to-noise ratio across the whole wavelength 
range from 3650 to 8500 A. Also shown for reference in Fig.|S| 
are the central band-passes of the 23 spectral indices of the 
Lick system (blue bands). Several of these, such as H<5a, 
Ca4227, G4300, H7A, H/3, Ca4455 and NaD, coincide with 
spectral features identified by MOPED to be most sensitive 
to the fractions of stars in various age ranges. A full coin- 
cidence is not expected between the red histogram and the 
blue bands, because many Lick indices are expec ted to be 
more sensitive to metallicity rather than age fe.g.. lWorthevl 
Il994h . 

Several spectral regions identified by MOPED to be 
highly sensitive to age in Fig. |S| do not coincide with any 
Lick index. To highlight all the features expected to reveal 
the presence of stars in various age ranges, we smooth each 
red histogram of Fig. [5] using a Gaussian kernel of rms 10 A 
and show the result as a solid curve in Fig. |7] For each age 
range, we mark with blue bands the spectral regions where 
the components of the smoothed vector exceed 1.5 times the 
rms variation over the 3650-8500 A wavelength range (we 
choose this level empirically to highlight only a moderate 
number of regions). Most of these regions appear to lie in 
the blue part of the optical spectrum; they could constitute 
the central band-passes of a set of optimal indicators of the 
presence of stars in various age ranges, which could be ex- 
ploited without appealing to the full MOPED analysis. In 
practice, only a few of the most sensitive indices should be 
selected to trace the presence of stars in each age range. We 
plan to explore this and the dependence of the spectral indi- 
cators on velocity dispersion and metallicity in more detail 
in a forthcoming analysis. 
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Figure 5. Relative contributions by stars in different age bins to 
the integrated spectra of 13.5 Gyr-old model galaxies, for 3 differ- 
ent time scales of star formation, r = 2, 7 Gyr and 00 (constant 
star formation rate), and for solar metallicity. For r = 7 Gyr, stars 
with ages 0.5-4 Gyr contribute only little to the integrated light at 
any wavelength. This makes the mass fraction of intermediate-age 
stars difficult to constrain on the basis of the integrated spectra of 
galaxies with smoothly declining star formation rates (see Fig.l2"l. 



4 APPLICATION TO THE INTERPRETATION 
OF OBSERVED GALAXY SPECTRA 

In this section, we show examples of derived star formation 
histories from observed, medium-resolution galaxy spectra 
using the method outlined in Section [5] We first recall the 
main characteristics of the SDSS EDR spectra we consider. 
We then show examples of derived star formation histories 
of early-type and late-type galaxies of comparable mass. Fi- 
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Figure 6. Components of the MOPED weighting vectors (red histogram) quantifying the sensitivity of the flux at each 
wavelength to the mass fractions of stars in the 6 age bins used in the present analysis, as indicated in the legend. For 
illustrative purposes, the amplitudes have been renormalized in each panel. Also shown for reference are the central band- 
passes of the 23 spectral indices of the Lick system (blue bands). The predictions are for a stellar velocity dispersion 
ay = 100 km s — 1 , solar metallicity and a fixed signal-to-noise ratio per pixel across the whole wavelength range from 3650 
to 8500 A. 
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Figure 7. Components of the MOPED weighting vectors quantifying the sensitivity of the flux at each wavelength to 
the mass fractions of stars in the 6 age bins used in the present analysis (as indicated), after smoothing with a Gaussian 
kernel of rms 10 A (solid curve). For illustrative purposes, the amplitudes have been renormalized in each panel. The blue 
bands mark spectral regions where the components of the smoothed vector exceed 1.5 times the rms variation over the 
3650—8500 A wavelength range. The prediction are for a stellar velocity dispersion ov = 100 km s — 1 , solar metallicity and 
a fixed signal-to-noise ratio across the whole wavelength range from 3650 to 8500 A. 
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nally, we constrain the global assembly of stellar mass in 
galaxies in different mass ranges in the Universe. 

4.1 SDSS galaxy spectra 

The observati onal sample we consid er is the SDSS Early 
Data Release JStoughton et al.ll2002T) . The SDSS is obtain- 
ing u, g, r, i, and z photometry of almost a quarter of 
the sky and spectra of at least 700,000 objects. The 'main 
galaxy sample' of the EDR includes the spectra of 32,269 
galaxies with r-band Petrosian magnitudes in the range 
14.50< r < 17.77, after correc tion for foreground Galactic ex- 
tinction iStrauss et aljl2002f) . The median redshift is about 
0.1. The SDSS spectra are acquired using 3-arcsecond di- 
ameter fibres that are positioned as close as possible to 
the centres of the target galaxies. The spectra are flux- 
and wavelength-calibrated, with 4096 pixels from 3800 A 
to 9200 A at resolving power A/ A A £s 1800. This is simi- 
lar to the resolution of our models in the wavelength range 
from 3650 A to 8500 A. We stress that, since we are inter- 
ested in only the high-pass signal of the spectra, our results 
are not sensitive to uncertainties in the spectro-photometric 
calibration of the SDSS EDR spectra ( Section YT%. 

We consider here the spectra of all galaxies with red- 
shifts in the range 0.005 < z < 0.2. We use smear- 
correcte d observations reduce d through the spectro-ld 
pipeline iSubbaRao et al J2002D . which we convert from vac- 
uum to air for comparison with the models. We remove 
sky (emission and absorption) lines. We also remove all 
detectable nebular emission line s using the procedure de- 
scribed in lTremonti et all 12004^1 that is based on accurate 
fits to the emission line-free regions of the spectra with 
model spectra broadened to the observed stellar velocity dis- 
persion. Finally, we mask those regions of the spectra that 
are known to strongly depend on changes in metal abun- 
dance ratios (Section l2.21 . 

4.2 The star formation histories of early-type and 
late-type galaxies 

To illustrate the constraints derived from high-pass spec- 
tra on the star formation histories of galaxies, we differen- 
tiate between early-type and late-type galaxies. We select 
the two types of galaxies on the basis of the 'concentration 
parameter', defined as the ratio C = R90/R50 of the radii 
enclosing 90 and 50 percent of the Petrosian r-band luminos- 
ity of a galaxy. This has been show n to correlate well with 
'by eye' morphological classification dShimasaku et aJl200li 
IStrateva et aljl200ll) . lstrateva et all l l200lh propose a cut at 
C — 2.6 to separate early- from late-type galaxies. We clas- 
sify conservatively as 'early type' all galaxies with C > 2.8 
and as 'late type' all galaxies with C < 2.4. We exclude 
gal axies identified as hosts of active galactic nuclei (AGNs) 
bv lKauffmann et alJ fcOOSt . These various cuts reduce our 
sample to 3229 early-type and 9430 late-type galaxies with 
redshifts 0.005 < z < 0.2. 

We derive star formation histories for all these galax- 
ies using the method described in Section [5] In this pro- 
cedure, we always adopt models with stellar velocity dis- 
persions close to the observed one and compute MOPED 
weighting vectors specific to each galaxy, based on the 



observed noise spectrum (Section 12.4.11 . As an example, 
we compare in Figs |H| and [§] the observed high-pass spec- 
tra of an early-type galaxy (#274-51913-392, IAU: SDSS 
J103840.25+002251.6) and a late-type galaxy (#268-51633- 
202, IAU: SDSS J095645. 76+011448.7) with the spectra of 
the models corresponding to the most likely star formation 
histories selected by our analysis. The early- type galaxy 
has C = 2.95, a 4000 A discontinuity D400 = 1.97 (as 
comp uted using the narrow-band definition of Balogh et al. 
1999), a velocity dispersion av = 127km s _1 and a me- 
dian signal-to-noise ratio per pixel S/N « 27. The late- type 
galaxy has C = 2.39, D4000 = 1.13, cry < 75 km s" 1 and 
S/N « 22. There is good agreement between models and 
observations for the main stellar (Balmer and metallic) ab- 
sorption features in the spectra of both galaxies. 

Figs 1101 and 1111 show star formation histories derived 
from the high-pass spectra of the early-type and late-type 
galaxies in Figs |H| and [5] respectively. In each case, the left- 
hand panel shows the evolution of the star formation rate, 
while the right-hand panel shows the cumulative assembly of 
stellar mass. The star formation rates in Figs llOl and llll were 
computed by dividing the mass fraction of stars formed in 
each bin of look-back time by the duration of that bin. Thus, 
they are normalized to a total of 1 Mq of stars formed over 
all ages. As in Figs|5J2| the thick line, box and error bars 
for each bin of look-back time in Figs 1121 and 1131 indicate 
the median value, 16-84 and 2.5-97.5 percentile ranges of 
the quantity plotted on the y-axis (these ranges indicate the 
68 and 95 percent confidence intervals in the recovered star 
formation history; the results of Section ^. 2l suggest that the 
16-84 percentile range is typical of the accuracy to which our 
method can recover star formation histories from high-pass 
galaxy spectra). Fig. 1101 shows that most of the stars in the 
early-type galaxy SDSS J103840.25+002251.6 are likely to 
have formed over 8 Gyr ago, while star formation appears 
to have occurred much later in the late-type galaxy SDSS 
J095645. 76+011448. 7. The significant recent star formation 
in this galaxy makes the mass fraction of stars with ages 
around 1 Gyr difficult to constrain (see also Section f3. 31 . 

It is interesting to compare the typical star formation 
histor ies of early-type and late - type galaxies of similar stellar 
mass . IKauffmann et al.l J2003T) derived stellar masses for the 
galaxies in the SDSS EDR, based on the observed strengths 
of D4000 and of the H<5a absorption line. We follow their 
procedure here and compute, for each galaxy in our sam- 
ple, the median mass-to- light ratio in the observer- frame i 
band that is predicted by our analysis of the high-pass spec- 
trum. 6 We assume that this ratio is the same for the whole 
galaxy as it is in the region sampled by the fibre. We then 
compute the total stellar mass A/* from the observed i-band 
luminosity, after correcting the i-band flux for attenuation 
by dust and the contaminatio n by emission lines, usin g the 
same procedure as outlined in Kauffma nn et alJ (|2003j). The 
galaxy stellar masses obtained in this way agree to withi n 
50 percent with those derived bv iKauffmann et al.l (2003). 



6 IKauffmann et ail i2003T) used the z band in their analysis, 
which falls outside the spectral range of the medium-resolution 
models considered here for galaxies at redshifts z < 0.2 (Sec- 
tion ooj. 
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Figure 8. Comparison of the high-pass spectrum of the early-type SDSS galaxy #274-51913-392 (IAU: SDSS 
J103840. 25+002251. 6, in blue) with that of the model corresponding to the most likely star formation history selected 
by the method described in Section [3] (in red). 




6500 



7000 



7500 6000 
Wavelength [ A 1 



8500 



9000 



Figure 9. Comparison of the high-pass spectrum of the late-type SDSS galaxy #268-51633-202 (IAU: SDSS 
J095645. 76+011448. 7, in blue) with that of the model corresponding to the most likely star formation history selected 
by the method described in Section [5] (in red). 



For the SDSS-EDR sample, the distribution of galaxy stellar 
masses peaks around 3 x 10 10 Mq . 

We focus on the typical star formation histories derived 
from the high-pass spectra of early-type and late- type galax- 
ies around the peak of the stellar mass distribution, i.e. with 
masses in the range 9.5 < log(M*/M0) < 11.5. One way to 
constrain these is to stack the spectra of all galaxies of each 
type in this mass range and to derive the corresponding star 
formation history from the stacked spectrum. In practice, 
we first arrange the galaxies in 4 bins of velocity dispersion 



(cry), 5 bins of 4000-A discontinuity (D4000) and 2 bins of 
redshift (z). Then, within each 3-dimensional bin, we stack 
the (emission-line subtracted) spectra of all galaxies. In this 
procedure, we normalize each spectrum to the total absolute 
magnitude derived from the observed Petrosian r-band mag- 
nitude. We also weight each spectrum by a factor 1/VJnax, 
where Kn ax is the largest volume within which the galaxy 
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Figure 10. Evolution of the star formation rate (left; normalized 
to 1 Mq of stars formed over all ages) and cumulative assembly 
of stellar mass (right) derived from the high-pass spectrum of the 
galaxy in Fig. [S] For each bin of look-back time, the thick line, 
box and error bars indicate the median value, 16-84 and 2.5- 
97.5 percentile ranges of the quantity plotted on the y-axis (these 
ranges indicate the 68 and 95 percent confidence intervals in the 
recovered star formation history). Bins with no star formation 
are shown blank, as is the last bin of the cumulative assembly of 
stellar mass (which is unity by construction). 
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Figure 11. Evolution of the star formation rate (left; normalized 
to 1 Mq of stars formed over all ages) and cumulative assembly 
of stellar mass (right) derived from the high-pass spectrum of the 
galaxy in Fig. [5] For each bin of look-back time, the thick line, 
box and error bars indicate the median value, 16-84 and 2.5- 
97.5 percentile ranges of the quantity plotted on the j/-axis (these 
ranges indicate the 68 and 95 percent confidence intervals in the 
recovered star formation history). Bins with no star formation are 
shown blank. 



would pass the sample selection criteria. 7 We derive the star 
formation history from the high-pass signal in the stacked 
spectrum in each bin of ov, -D4000 and z. This star for- 
mation history corresponds to a final stellar mass density of 
£i(M*,j/Vmax,i). The typical star formation history of galax- 
ies of a given type is therefore given by the Ej(M,,i/V^ naXj i)- 
weighted average of the star formation histories of all bins 
for that type. 

Figs ll2l and ll3l show the typical star formation histories 
derived in this way from the high-pass spectra of 659 early- 
type and 385 late-type galaxies with stellar masses in the 
range 9.5 < Iog(M*/Mo) < 11.5 (i.e., around the peak of 
the stellar mass distribution of SDSS-EDR galaxies) in our 
sample. In each case, the left-hand panel shows the typical, 
class-averaged evolution of the star formation rate, while 
the right-hand panel shows the cumulative assembly of stel- 
lar mass. Figs ll2l shows that most of the stars in the typical 
early-type galaxy in our sample are likely to have formed 
over 8 Gyr ago, although a small fraction of the total stel- 
lar mass can be accounted for by stars with ages down to 
4 Gyr. This is consistent with a wide range of constraints on 
t he fraction of intermediate-age stars in early-type galaxies 
Jj0rgensenlll999l:lTrager et aLlboOfj iMenanteau et alJl200lt 
iBell et al.ll2003l) . We note that iReichardt et alJ <200ll) used 
MOPED to interpret the low-resolution spectra of 8 bright 
nearby E/SO galaxies from the Kennicutt (1992) atlas. They 
found that stars younger than 3 Gyr could account for over 
30 percent of the total stellar mass for 6 of these galaxies. 
This is not the case for the typical early-type galaxy in our 
sample. 

Fig. 1131 shows how the typical star formation history of 
a late-type galaxy differs from that of an early-type galaxy 
in the same stellar-mass range. The typical late- type galaxy 
with 9.5 < log(M*/M0) < 11.5 appears to have formed 
stars at a roughly constant or smoothly declining rate over 
the lifetime of the universe. The apparent 'drop' in the me- 

7 We compute V max using the effective area and the bright (14.50 
mag) and low (17.77 mag) magnitude cutoffs of the SDSS EDR. 
We do not include any evolutionary corrections to Vmax. 



dian star formation rate at look-back times around 1 Gyr 
in Fig 113b is not statistically significant. It is a likely con- 
sequence of the difficulty of constraining the mass fraction 
of stars of this age for smoothly declining star formation 
histories (Section 13.31 . We note that the average evolution 
of the star formation rate of a late-type galaxy in Fig. 113b 
is consistent wi th the star formation his tory of the Milky 
Way derived bv lRocha-Pinto et al.l i2000l) from the analysis 
of 552 Galactic stars. 

The results of Figs 1121 and 1131 illustrate the accuracy 
to which our method can recover star formation histories 
from the observed high-pass spectra of individual galaxies 
in the SDSS EDR. The typical uncertainties in the derived 
star formation histories are somewhat larger than those ex- 
pected from the analysis of model galaxy spectra with simi- 
lar signal-to-noise ratio (S/N = 30) in Sectional (FigsUgJ. 
This arises in part from our simplifying assumption that all 
stars in the model galaxies studied in Section ^. 2l had a fixed 
metallicity (corresponding to the effective metallicity). The 
mixing of stars with different chemical compositions in real 
galaxies makes the interpretation of the high-pass spectrum 
using single-metallicity models less precise. Also, our mod- 
els have fixed, solar metal abundance ratios, whereas the 
abundance ratios of heavy elements are expected to vary in 
external galaxies (Section 12. 21 . The influence of this uncer- 
tainty on the present study should be minimized by the fact 
that we have masked those features that are expected to be 
most sensitive to changes in metal-abundance ratios in the 
SDSS galaxy spectra. 

We have also derived rough constraints on the metal- 
licities of all galaxies in our sample by assigning them sub- 
solar {QAZq), solar (Zq) or supra-solar (1.5-Zq) metallicity, 
as described in Section ^. 4. 21 We have found that the results 
compare relatively well t o the more refined metallicities de- 
rived bv lGallazzi et all (I2004T) for these galaxies, based on 
the analysis of selected Balmer and metallic lines with weak 
dependence on metal-abundance ratios. 
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Figure 12. Typical star formation history (left; normalized to 
1 Mq of stars formed over all ages) and cumulative assembly of 
stellar mass (right) derived from the high-pass spectra of early- 
type galaxies with stellar mass in the range 9.5 < log(M,/MQ) < 
11.5 in the SDSS EDR. The results are based on the spectra of 
659 galaxies with concentration parameters C > 2.8 (see text for 
detail). The boxes have the same meaning as in Fig. 1101 
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Figure 13. Typical star formation history (left; normalized to 
1 Mq of stars formed over all ages) and cumulative assembly of 
stellar mass (right) derived from the high-pass spectra of latc- 
type galaxies with stellar mass in the range 9.5 < log(M„ /Mq) < 
11.5 in the SDSS EDR. The results are based on the spectra of 
385 galaxies with concentration parameters C < 2.4 (see text for 
detail). The boxes have the same meaning as in Fig. I1UI 



4.3 Cosmic star formation history 

We now use the high-pass spectra of SDSS galaxies to set 
constraints on the global star formation history of the Uni- 
verse and its distribution among galaxies in different mass 
ranges. Here, we include all galaxies at redshifts 0.005 < z < 
0.2 in the SDSS EDR sample. We also include those g alax- 
ies identified as AGN hosts bv lKauffmann et al.l i2003t) . As 
shown by these authors, the presence of a (type 2) AGN does 
not alter significantly the stellar absorption-line signatures 
in the host-galaxy spectrum. By analogy with our approach 
in Section 14.21 we proceed by first stacking the spectra of 
galaxies in different bins of velocity dispersion (ov; 4 bins), 
4000-A discontinuity (D4000; 3 bins) and redshift (z; 2 bins) 
to improve the signal-to-noise ratio. As before, we normal- 
ize each individual spectrum to the total absolute magni- 
tude derived from the observed Petrosian r-band magnitude 
of the galaxy and apply a weight of 1/Knax- Then, we de- 
rive the star formation history from the high-pass signal in 
the stacked spectrum in each bin of ay, 734000 and z. This 
yields the probability density functions of the mass fractions 
of stars formed in 6 age bins for this spectrum. This star for- 
mation history corresponds to a final stellar mass density of 
E^M^i/Knaxj), which we use to transform stellar mass frac- 
tions into stellar mass densities. To obtain the global star 



formation history per unit comoving volume of the Universe, 
we then simply combine the star formation histories over all 
bins of ov, -D4000 and z. In practice, this is obtained by com- 
bining the probability density functions of the stellar mass 
densities formed in each r edshift bin, following the precepts 
outlined in appendix A of iBrinchmann et alJ l!20o3) . 

Fig. 1141 shows the evolution of the star formation den- 
sity (left) and cumulative assembly of stellar mass (right) 
as accounted for by galaxies in different mass ranges in a 
unit comoving volume of the Universe: (a) 10 7 -10 9 Mq ; (b) 
1O 9 -1O U M ; (d) 10 11 - 10 13 M Q ; and (d) full mass range 
(10 7 -10 13 Mq). The most striking result from this figure is 
the fact that the stars that end up in massive galaxies today 
appear to have formed on average earlier than those that end 
up in low-mass galaxies. This result was alread y obtained 
in a different analysis by iHeavens et alj i2004l) , who used 
MOPED to interpret the low-resolution spectra of 96,545 
galaxies in the SDSS Data Release One (DR1). It is consis- 
tent with the 'down-sizing' scenario, in which the mass and 
luminosity of the galaxies undergoing active star formation 
become progressive ly lower as the Universe becomes older 
iCowie et alJll996[). We note that bo th the present analy- 
sis and that of IHeavens et ail i2004l) rely on the spectral 
analysis of nearby galaxies. Hence neither analysis indicates 
whether the stars formed once the galaxies were already as- 
sembled into systems with masses similar to the present-day 
ones. 

It is unlikely that the results of Fig. 1141 be seriously 
affected b y the ma g nitud e limit of our sample. As shown by 
IBrinchmann et al. ( 2004}) , at the low-redshift limit of z = 
0.005, the r = 17.77 magnitude limit of the spectroscopic 
SDSS sample corresponds roughly to a stellar mass limit of 
just under 1O 8 M0 for a galaxy which would have formed all 
its stars 13.5 Gyr ago and faded ever since. Most galaxies of 
the same mass containing younger stars would be brighter, 
as attenuation by dust is not ex pected to be significant i n 
low-mass galaxies (see figure 6 of IBrinchmann et al.ll2004l) . 
We find that restricting the stellar-mass range of galaxies in 
Fig. 114b to the more conservative completeness range 8 < 
log(M*/M0) < 9 has a negligible effect on the derived star 
formation history. Thus, the magnitude limit of the SDSS 
sample cannot account for the conclusion from Fig. I14l that 
the stars in massive galaxies today formed on average earlier 
than those in low-mass galaxies. 

Aperture effects could potentially affect the results of 
Fig. ll4l because. at the low redshifts of the SDSS sample, the 
small-aperture SDSS fibres tend to sample the inner, bulge- 
dominated re^ioirs_ofjnassive_galaxies wit h bulge+disc mor- 
pholo gies. iKewlev et al.l i2005l see also iGlazebrook et al] 
2003) find that the influence of aperture effects on star- 
formation-rate estimates is significantly reduced when at 
least 20 percent of a galaxy's light is gathered in the fibre. 
They show that, for the average SDSS galaxy, this corre- 
sponds to restricting studies to redshifts larger than z ~ 
0.04. Since most galaxies more massive than M* ~ 1O 9 M0 
in our sample lie at redshifts greater than this limit, aper- 
ture effects should have only a weak influence on the results 
of Figs ll4b and c. Low-mass galaxies in Fig. 114b are not ex- 
pected to be subject to strong aperture effects either, since 
their morphologies tend to be dominated by either a bulge 
or a disc, but not both. T his is supported by the finding by 
IBrinchmann et all J2004I) that low-mass SDSS galaxies show 
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Figure 14. Evolution of the star formation density (left) and 
cumulative assembly of stellar mass (right) derived from the high- 
pass spectra of SDSS galaxies in different mass ranges at redshifts 
0.005 < z < 0.2. (a) For galaxies with stellar masses in the range 
7 < log(M*/M0) < 9. (b) For galaxies with stellar masses in 
the range 9 < log(M*/MQ) < 11. (c) For galaxies with stellar 
masses in the range 11 < log(M*/M@) < 13. (d) For galaxies 
with stellar masses in the full range 7 < log(M t /Mq) < 13. For 
each bin of redshift, the thick line, box and error bars indicate 
the median value, 16-84 and 2.5—97.5 percentile ranges of the 
quantity plotted on the y-axis (these ranges indicate the 68 and 
95 percent confidence intervals in the recovered star formation 
history). Bins with no star formation are shown blank, and the 
last bin of the cumulative assembly of stars is purposely omitted 
as it is unity by construction. 



only weak colour gradients (we note that the star formation 
density in the lowest redshift bin in Fig. ll4H . which is based 
on the spectra of only ~ 10 nearby low-mass galaxies, is 
highly uncertain). Thus, we conclude that aperture effects 
should not strongly influence the results of Fig. 1141 

The evolution of the total star formation density 
of the Univer s e in Fig. I14H differs from that found by 
Irleavens et alJ <l2004t) in that we do not find any statisti- 
cally significant peak in the star formation activity at red- 
shift 2 <; 1. This results in part from our more conserva- 
tive treatment of the uncertainties associated to the star 
formation histories derived with MOPED (see Section [2.411 . 
The relatively large uncertainties in Fig. 1141 also reflect the 
limitations of current models for interpreting the high-pass 
spectra of observed galaxies (see the discussion at the end 
of Section f4.21 . These uncertainties should be taken as very 
conservative, in the sense that they were constructed ig- 
noring any possible correlation between the star formation 
histories of galaxies in different bins of velocity dispersion, 
4000-A discontinuity and redshift classes (see above). Ac- 
cording to Fig. I14H . about 40 ± 5 percent of all stars in 
the Universe formed prior to redshift z » 1, and the star 
formation density has since been declining sm oothly. These 
result s are consistent with those obtained bv lBaldrv et alJ 
(2002) from the analysis of the stacked, low- resolution spec- 
tra of 166,000 galaxies at redshifts 0.03 < z < 0.25 ex- 
tracted fr om the Two-degree Field Galaxy Redshift Survey 
(2dFGRS: IColless et al1l200ll) . They conclude that the star 
formation density of the Universe must have peaked in the 
past (z > 0.6) at a value at least three times the present-day 
one, while no more than 80 percent of all stars must have 
formed at redshifts z > 1. 



5 SUMMARY AND CONCLUSIONS 

We have presented a method for extracting star formation 
histories from medium-resolution galaxy spectra, which is 
free of uncertainties arising from the spectro-photometric 
calibration and intrinsic attenuation by dust in the galax- 
ies. The method focuses on the interpretation of the high- 
pass signal in galaxy spectra, using a combination of the 
data compression algorithm MOPED jHeavens et al]l2000ft 
with the medium-r e soluti on population synthesis code of 
iBruzual fc Charlod l|2003l) . which covers the wavelength 
range 3650-8500 A at a resolving power of ~ 2000. The 
principle of this method is to quantify the sensitivity of each 
pixel in a high-pass galaxy spectrum to the mass fraction of 
stars formed in 6 age bins spanning the age of the Uni- 
verse. Given an observed spectrum, statistical estimates of 
the mass fraction of stars formed in each age bin allow one 
to constrain the entire star formation history of the galaxy. 
Our implementation of MOPED differs in several ways from 
previous applications of this algorit hm to the interpreta- 
tion of low-resoluti on galaxy spectra l)Reichardt et al.ll200ll : 
iPanter et alJ I2OO3T) . Overall, these differences should lead 
to more conservative estimates of the errors in the derived 
star formation histories (Section 12.41 . Our approach does 
not require prior knowledge of the metallicity but assumes 
that all stars in a given galaxy have the same metallicity, 
which should be regarded as the effective metallicity (i.e., 
the metallicity of the stars dominating the light). 
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We have tested in detail the ability of our method to 
recover star formation histories from model high-pass spec- 
tra of galaxies with different time-scales of star formation. 
We find that the method can recover the full star forma- 
tion histories of these models, provided that the median 
signal-to- noise ratio per pixel is large, i.e. S/N > 20-30. 
Our method allows us to identify those spectral features 
that are most sensitive to the presence of stars of a given 
age in the integrated spectrum of a galaxy. We have ex- 
emplified how this could be used to define new 'Lick-type' 
spectral indices for studies of the star formation histories 
of galaxies (Section 13.41 . The application of our method to 
model galaxy spectra also reveals a fundamental limitation 
in the recovery of the star formation histories of galaxies for 
which the optical signatures of intermediate-age stars are 
masked by those of younger and older stars. One of the con- 
sequences of this limitation is that the mass fraction of stars 
with ages around 1 Gyr is difficult to recover with accuracy 
in galaxies with smoothly declining star formation histories 
(i.e. with star formation time-scales of about 7 Gyr). Inter- 
estingly, this type of star formation history i s expected to 
be re presentative of the Universe as a whole fealdrv et alJ 
120021) . 

As an illustration, we have used our method to derive 
star formation histories from high-quality galaxy spectra in 
the SDSS EDR. To perform this analysis, we have first re- 
moved emission lin es from the observed s pectra using the 
method outlined bv lTremonti et ail (I2004T) . We have found 
that the star formation histories of morphologically identi- 
fied early- and late-type galaxies with stellar masses near 
the peak of the distribution in the local Universe, i.e. with 
9.5 < log(M*/M0) < 11.5, are in agreement with what is 
commonly assumed about these systems. Early-type galax- 
ies appear to have formed most of their stars over 8 Gyr ago, 
although a small fraction of the total stellar mass of these 
galaxies may be accounted for by stars with ages down to 
4 Gyr. In contrast, late- type galaxies appear to have formed 
stars at a roughly constant rate. We also examined the con- 
straints set by the high-pass signal in the stacked spectra of 
a complete magnitude-limited sample of 20,623 SDSS-EDR 
galaxies on the global star formation history of the Universe 
and its distribution among galaxies in different mass ranges. 
Our results confirm that the stellar populations in the most 
massive galaxies today appear to have for med on average 
earlie r than those in the least massive ones jHeavens et alJ 
120041) . Our analysis does not reveal any statistically signifi- 
cant peak in the star formation activity of the Universe at 
redshifts below unity. 

At present, the main limitation in the interpretation of 
medium-resolution galaxy spectra is the fixed, solar metal 
abundance ratios of the models at all metallicities while the 
abundance ratios of h eavy elements are observed to vary in 
exter nal galaxies fe-g- JWorthev et al.ll 992: Eiscnstcin et alJ 
2003). To minimize the consequences of this limitation, we 
have masked those regions of galaxy spectra that are known 
to strongly depend to changes in metal abundance ratios 
(Section 12. 21 . However, we suspect that at least part of 
the relatively large uncertainties affecting the star formation 
histories derived from the high-pass spectra of SDSS-EDR 
galaxies in Section 2] may result from this limitation. On 
the theoretical front, work is underway to develop the abil- 
ity to model the spectral evolution of galaxies with various 



metal abundances ratios at medium resolution (e.g., Coelho 
et al., in preparation). It is worth noting that the (weak) 
constraints derived from our analysis on the metallicities of 
SDSS-EDR galaxies compare relatively well t o the more re- 
fined ones derived for the same galaxies by iGallazzi et alJ 
i2004h from the analysis of selected Balmer and metallic 
lines with weak dependence on metal-abundance ratios (Sec- 
tion l4"2l. 

The results obtained in this paper open the door to fur- 
ther refined interpretations of galaxy spectra in terms of star 
formation history, metallicity and dust content. For exam- 
ple, the low-pass signal in galaxy spectra, which is deliber- 
ately left aside in our analysis, contains valuable information 
about intrinsic attenuation by dust. In Section ^. 3. 31 we de- 
scribed how this information can be extracted from the low- 
pass signal, based on the constraints set by the high-pass 
signal on the star formation history (this requires that the 
spectra be spectro- photometrically well calibrated). Another 
interesting development should come from the application of 
our method to the interpretation of galaxy spectra gathered 
by high-redshi ft surveys, suc h as the VIRMOS-VLT Deep 
Survey <ILe Fevre et al j |2001f) an d the Deep Extragalactic 
Evolutionary Probe IIDavis et all l2003L Such analyses will 
allow a more precise accounting of the star formation history 
of galaxies at the epoch z > 1, which is poorly sampled by 
the analysis of local galaxies. These various developments 
should lead us to a deeper understanding of how galaxies 
formed and evolved. 
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